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i Abstract 

a.: 

r-| ! The photoproduction of prompt photons, together with an accompanying jet, has 

been studied in ep collisions at a centre-of-mass energy of 318 GeV with the ZEUS 
detector at HERA using an integrated luminosity of 77 pb" 1 . Cross sections were 
measured for the transverse energy of the photon and the jet larger than 5 and 
6 GeV, respectively. The differential 7+jet cross sections were reconstructed as 
functions of the transverse energy, pseudorapidity and x° bs , the fraction of the 
incoming photon momentum taken by the photon-jet system. Predictions based 
on leading-logarithm parton-shower Monte Carlo models and next-to-leading- 
order (NLO) QCD generally underestimate the cross sections for the transverse 
energies of prompt photons below 7 GeV, while the fc^-factorisation QCD cal- 
culation agrees with the data better. When the minimum transverse energy of 
prompt photons is increased to 7 GeV, both NLO QCD and the /cT-factorisation 
calculations are in good agreement with the data. 
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1 Introduction 



Events containing an isolated photon (prompt photon) are a powerful tool to study hard 
interaction processes since such photons emerge without the hadronisation phase by which 
a final state quark or gluon forms a jet. In ep collisions, the presence of a jet in addition 
to the photon allows measurements that are more sensitive to the underlying partonic 
processes than is possible for inclusive prompt-photon events. In particular, final states 
with a prompt photon with a high transverse energy (Et) together with a high-ET jet 
are directly sensitive to the quark content of the proton through the scattering of the 
exchanged photon with a quark, 75 — > 75 (Compton scattering). In this case, the inci- 
dent photon is point like, and the process (Fig. [T] (a-b)) is called direct. For exchanged 
four-momentum transfer close to zero (photoproduction), the additional contribution to 
prompt-photon events from the gq — > q'y process, in which one of the initial partons comes 
from the photon which displays a hadronic structure (resolved process, Figs. U(c-d)), can 
be dominant [1-3]. Prompt-photon measurements can be used to constrain the parton 
distribution functions (PDFs) of the proton and of the photon, as well as provide a testing 
ground for QCD calculations. A number of predictions exist [1-4] that can be confronted 
with the data. 

The first observation by ZEUS of isolated photons accompanied by a hadronic jet in 
photoproduction used an integrated luminosity of 6.4 pb _1 [5]. Distributions sensitive to 
the intrinsic kx in the 7+jet final state were later measured by ZEUS [6]. Inclusive prompt- 
photon cross sections with no jet requirement have also been reported [7]. Recently, HI 
have published results on the 7+jet final state in photoproduction [8]. 

This paper reports the first ZEUS results on differential cross sections of the 7+jet final 
state in the photoproduction regime of ep scattering. The cross sections are presented 
as a function of the transverse energy and pseudorapidity of both the photon (Ej,, rj 1 ) 
and the jet (E^ , if et ), as well as the fraction of the incoming photon momentum taken 
by the photon-jet system (x° bs ). In contrast to previous measurements [5-8], the present 
analysis is based on the conversion-probability method, which uses information on the 
frequency with which photons convert to e + e~ in front of a dedicated preshower detec- 
tor. Cross sections for 7+jet events are compared to next-to- leading-order (NLO) QCD, 
calculations based on the /cr-factorisation approach and Monte Carlo (MC) models incor- 
porating leading-order matrix elements plus parton showers. Since jets at relatively low 
E^ are measured in addition to the photon, parton-level calculations were corrected for 
hadronisation effects using a MC model. The hadronisation correction for 7+jet events is 
expected to be smaller than for dijets with similar jet transverse energies since the photon 
does not undergo hadronisation. Therefore, for low E^ , the theoretical predictions for 
the 7+jet cross sections are expected to be more reliable than for dijet final states. 
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2 Data sample and experimental setup 



The data sample was taken during the 1999-2000 period, corresponding to an integrated 
luminosity of 77.1 ± 1.6 pb" 1 . The positron or electron beam energy was 27.5 GeV and 
the proton beam energy was 920 GeV, corresponding to a centre-of-mass energy of 318 
GeV. Here and in the following, the term "electron" denotes generically both the electron 
(e - ) and the positron (e + ), unless otherwise stated. 

ZEUS is a multipurpose detector described in detail elsewhere [9]. Of particular im- 
portance in the present study are the central tracking detector, the uranium-scintillator 
calorimeter and the barrel preshower detector. 

The central tracking detector (CTD) [10] is a cylindrical drift chamber with nine super- 
layers covering the polar-angle 1 region 15° < 9 < 164° and the radial range 18.2 — 79.4 cm. 
Each super-layer consists of eight sense-wire layers. The transverse-momentum resolution 
for charged tracks traversing all CTD layers is cr{p T )/p T — 0.0058pr© 0.0065 ©0.0014/pr, 
with pt in GeV. 

The CTD is surrounded by the uranium-scintillator calorimeter, CAL [11], which is di- 
vided into three parts: forward, barrel and rear with the barrel consisting of 32 modules. 
The calorimeter is longitudinally segmented into electromagnetic (EMC) and hadronic 
(HAC) sections. The smallest subdivision of the CAL is called a cell. The energy resolu- 
tion of the calorimeter under test-beam conditions is oe/E = 0.18/ y/E for electrons and 
(Je/E = 0.35 /\[E for hadrons, with E in GeV. 

The luminosity was measured using the bremsstrahlung process ep — > epq with the 
luminosity monitor [12], a lead-scintillator calorimeter placed in the HERA tunnel at 
Z = -107 m. 

The ZEUS barrel preshower detector (BPRE) [13] is located in front of the barrel calorime- 
ter. The BPRE detector consists of 32 cassettes each containing 13 scintillator tiles of 
size 20 x 20 cm that were installed directly in front of each of the 32 barrel CAL modules. 
The measured output, calibrated in minimum ionising particle units (mips), is propor- 
tional to the energy loss of the incident particle after interaction with material (mainly 
the superconducting coil) in front of the barrel calorimeter. 

The mip calibration of each of the 416 channels of the BPRE was done using all events 
triggered in the ZEUS detector. A luminosity of approximately 1 pb _1 was required 
for each calibration run. The one-mip signal was validated using cosmic-ray muon data. 
The single-mip resolution was measured to be 0.3 mips and a minimum charge threshold 

1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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corresponding to this value was applied to each channel. After calibration and correction 
for dead or inefficient channels, the signal efficiency for scattered electrons from deep- 
inelastic events was larger than 99%. 

3 Theoretical predictions 

The measured 7+jet cross sections were compared to NLO QCD based on collinear factori- 
sation and DGLAP evolution [14], as well as to calculations based on the /^-factorisation 
approach with unintegrated quark and gluon densities. 

A NLO calculation with additional higher-order terms was performed by Krawczyk and 
Zembrzuski (KZ) [3]. The calculation includes the leading-order term jq — > jq, as 
corrections to this term, initial and final resolved-photon contributions, double-resolved 
contributions and the direct box diagram 7$ — > 7*7. The latter two contributions are 
calculated to order a 2 s . No intrinsic transverse momentum of the initial-state partons in 
the proton was assumed. The renormalisation and factorisation scales for such calculation 
are set to hr = \ip = Ej,. The GRV parameterisation of the proton PDF [15], the photon 
PDF [16] and the fragmentation function [17], were used. 

A similar NLO calculation by Fontannaz, Guillet and Heinrich (FGH) [2] contains addi- 
tional higher-order corrections to the resolved photon process. For the FGH calculation, 
the MRST01 [18] proton PDF and the AFG02 [19] photon PDF were used. 

The prediction of A. Lipatov and N. Zotov (LZ) [4] is based on the fey-factorisation [20] 
method. The LZ calculation uses the unintegrated quark and gluon densities of the 
proton and photon according to the Kimber-Martin-Ryskin (KMR) prescription [21] with 
the GRV parametrisations [15,16] of collinear quark and gluon densities. In this approach, 
both direct and resolved contributions are taken into account. 

For all the calculations discussed, jets were reconstructed by running the longitudinally 
invariant kp cluster algorithm in the inclusive mode [22] on partons. A prompt-photon 
jet was defined as a jet containing the final-state photon. An isolation requirement, 
^7, (true) ^ q g ^gj-g _£^>( true ) [ s the transverse energy of the final-state photon and Ej, 
is the total transverse energy of the prompt-photon jet, was applied to avoid the effects of 
collinear photon emission from quarks and to match the analysis isolation requirement (see 
Section EJ). A comparison with NLO calculations based on a cone isolation requirement 
showed consistent results [23]. 

The calculations were corrected for hadronisation effects using the Pythia MC model dis- 
cussed in Section HI These corrections, which are negligible in the case of inclusive prompt 
photons, cannot be neglected when an accompanying jet is required. The hadronisation 
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correction factors were defined as C had = ex (hadrons) /a (partons), where a denotes the 
differential cross sections calculated at the hadron and parton levels of the MC model, 
respectively. For both the parton and hadron levels of the MC generated events, the 
prompt photon was defined as the hr jet with at least one photon and with the isolation 
requirement E% > 0.9 Ej,. The E? and 77 distributions at the parton level in the MC 
model have a different shape than in the NLO calculations, especially at low Ej,, where 
the NLO predictions rise faster than do those of the MC. To determine the hadronisation 
corrections, the MC parton distributions were reweighted to match the shapes of the NLO 
calculations. The reweighting was performed in four dimensions defined by the Ej,, ?7 7 , 
E^ 1 and if et variables. 

The final hadronisation correction was determined from Pythia after the parton-level 
reweighting procedure discussed above. The Herwig model discussed in Section 0] re- 
quires a large reweighting so it was not used for the hadronisation correction. 

The hadronisation correction factor for the total cross section in the kinematic range 
defined in Section [7| was 0.92. The corrections for the 7 + jet differential cross sections 
are close to unity for large transverse momenta of the photon, but they decrease to 0.78 
at low Ej,. It was verified that if no jet was required, the hadronisation corrections were 
close to unity. 

The Pythia default setting includes a multiple-interaction simulation. It was verified that 
exclusion of multiple interactions from the parton-level of Pythia results in a negligible 
change in the hadronisation corrections. 

4 Monte Carlo simulation 

The measured cross sections were compared to leading-order Monte Carlo (MC) models 
which use the QCD parton shower approach to incorporate high-order QCD effects fol- 
lowed by fragmentation into hadrons. The MC events were generated with the Pythia 
6.3 [24-26] and with the Herwig 6.5 [27, 28] models using the default parameters in 
each case. The CTEQ5L [29] proton PDF was used together with the SaS-2D parame- 
terisation [30] for the photon PDF. Both direct and resolved prompt-photon events were 
generated. 

The same MC samples were used to calculate the acceptance and to evaluate the signal 
and background content of the sample. Samples of background photoproduction events 
(without prompt-photon subprocesses) were generated in addition to the prompt photon 
samples. Both direct and resolved processes were simulated. These MC samples provided 
background photons from the decay of hadrons (predominantly from 7r° mesons) . 
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The generated events were passed through a full simulation of the detector using Geant 
3.13 [31] and processed with the same reconstruction program as used for the data. The 
MC samples after the detector simulation do not give a good description of the Et and 
rj distributions seen in the data. Such discrepancies are most prominent at low E T , and 
were attributed to an inadequacy of the MC models. For the acceptance calculations, the 
MC distributions were reweighted to match the distributions in Et and rj of the data. The 
reweighting was performed in four-dimensional phase space in Et and 77 of the photon 
and of the accompanying jet; thus correlations between these kinematic variables were 
properly taken into account. 



5 Data selection and prompt-photon reconstruction 
5.1 Event selection and jet reconstruction 

The online selection made use of a standard ZEUS electron finding algorithm to select 
events with an electromagnetic cluster [5]. For the offline analysis, neutral-current deep 
inelastic (DIS) events with an identified scattered electron candidate were removed from 
the sample. This restricted the virtuality of the incident photon to the range Q 2 < 1 GeV 2 . 
In addition, the following cuts were applied: 

• charged current DIS events were rejected by requiring the total missing transverse 
momentum in the calorimeter to be less than 10 GeV; 

• 0.2 < i/jb < 0.8, where ?/jb is the inelasticity parameter reconstructed with the 
Jacquet-Blondel method [32]; 

• I Vertex |< 50 cm, where ^vertex is the event- vertex position determined from the 
tracks. 

Jets were reconstructed by running the longitudinally invariant kx algorithm in inclu- 
sive mode [22] on energy-flow objects (EFOs) [33], which are based on a combination of 
track and calorimeter information. The jet variables Et and 77 were defined according 
to the Snowmass convention [34]. Each jet was classified as either a photon candidate 
or a hadronic jet. The photon-candidate jet was required to consist of EFOs without 
associated tracks and to be within the CTD and BCAL acceptance, —0.74 < rp < 1.1. 
For this jet, E-EMc/E to t > 0.9 is required, where E^mc is the energy reconstructed in the 
electromagnetic part of the CAL and E tot is the total energy of this jet. After correction 
for energy losses, the cut E T > 5 GeV was applied. 

Hadronic jets, after correction for energy losses, were selected in the kinematic range 
&t > 6 GeV, -1.6 < rf ct < 2.4. They were required to have E^mc/E^ < 0.9. If more 
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than one jet was found within the above kinematic cuts, the jet with the highest Ej, was 
accepted. The minimum transverse momentum of the hadronic jet was set to be higher 
than for the photon candidate, since the NLO calculations employed in this analysis are 
unstable for symmetric cuts on the minimum transverse momenta of both jets [35,36]. 

In total, 3910 events with a prompt-photon candidate and a hadronic jet were selected. 

5.2 Identification of isolated photons and hadrons 

For the prompt-photon identification, the conversion-probability method based on the 
BPRE was used. In contrast to the shower-profile approach used in previous measure- 
ments [5-7], the present approach uses the probability of conversion of photons to e + e~ 
pairs in detector elements and inactive material, mainly the solenoid located in front of 
the BCAL. The conversion probability for a single photon is lower than for multiphoton 
events arising from neutral meson decays (tt°, rj, etc.); therefore, small BPRE signals can 
be used to identify isolated photons. 

The response of the BPRE to single isolated photons was studied using the deeply virtual 
Compton scattering (DVCS) data, ep — > e 7p, taken during 1999-2000. This sample is 
known to provide photons of high purity [37]. Events with two isolated electromagnetic 
clusters and one CTD track were pre-selected. One cluster was required to have energy 
E e i > 8 GeV and to be associated with the CTD track, thereby ensuring compatibility with 
the scattered electron. The cluster without an associated track was then reconstructed 
using the kr cluster algorithm as described in Section 15.11 The photon candidate was 
required to be in the BCAL region, —0.74 < rp < 1.1, and to have energy in the range 
5 < E 1 < 10 GeV. The BPRE signal for the DVCS photons was determined as the 
sum of the signal of the BPRE tiles whose centre falls within a cone of size 0.7 in r]-<j) 
around the photon candidate. A smaller cone size leads to an efficiency that is not well 
reproduced by the MC. Other details of the DVCS selection and MC simulations are given 
elsewhere [37,38]. 

In the DVCS sample, the fraction of events with BPRE signal below one mip is very sen- 
sitive to the amount of material in front of the BPRE since such events are dominated by 
non-converted photons. The MC simulation overestimates this fraction by 19% compared 
to the data due to an inadequate simulation of the material in front of the BPRE. This 
discrepancy does not have a significant dependence on the cone size. 

The amount of inactive material was further studied using scattered electrons from DIS 
events. This study indicated that more material in front of the BCAL was necessary for 
the MC simulation. 
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Using a dedicated Geant simulation [31], it was found that an increase of inactive ma- 
terial in the MC simulation by 0.25Xo was sufficient to describe the fraction of events 
without photon conversions seen in the DVCS data. For this Geant simulation, it was 
assumed that all inactive material is distributed uniformly in the region of the ZEUS 
solenoid located in front of the BPRE tiles. The effect of the additional material was 
then taken into account in the standard ZEUS MC by applying a correction to the BPRE 
distribution based on the results of the dedicated simulation. 

Figure El shows the BPRE signal for the DVCS data compared to the DVCS Monte Carlo 
model [38] after correction for additional material. There is good agreement between the 
data and the MC distribution. This shows that the inactive material and the BPRE 
resolution are well represented in the MC simulation. 

The MC response of the BPRE to single hadrons was studied for tt° and rj mesons. Since 
these mesons decay to several photons, more conversions to e + e~ will occur than for single 
photons. As expected, the average BPRE signal for 7r° and r\ mesons was larger than for 
the isolated photons. An example of the MC simulation of the BPRE response to single 
7, 7T° and r\ is shown in Figure El The BPRE distributions for ir° and r\ mesons were also 
corrected to take into account additional dead material in front of the BPRE. 

5.3 Extraction of the prompt-photon signal 

The BPRE signal for prompt-photon candidates selected as described in Sect. 15.11 was 
determined using a cone of radius 0.7 in r)-(f) space, as was done for the DVCS analysis. 
Figure 0] shows the comparison between the data and the Pythia MC for: the BPRE 
signal for the photon candidate; the difference between the total calorimeter energy and 
the energies of the jet and the photon candidate, AE = E tot — E^ ct — E 1 ; and the distance 
from the photon candidate to any EFO in an event, 



where r^ 7 (^efo) and 7 (0efo) are the azimuthal angle and pseudorapidity of the photon 
candidate (EFO). 

Figure El^a) shows that there is a significant fraction of events with a small number of 
mips, similar to the DVCS data. However, since the dijet photoproduction cross section 
is higher by several orders of magnitude than the 7+jet cross section, there is additional 
hadronic background even after the cuts discussed in Section 15.11 

The BPRE distribution for the prompt-photon candidates was used to determine the 
background fraction. The fraction of inclusive dijet photoproduction events needed was 
found from a ^-minimisation procedure. After the inclusion of the background events, 
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the shape of the BPRE distribution for the prompt-photon candidates is well reproduced 
by the MC simulation, as shown in Figure a). 

The inclusion of the dijet background leads to good description of the CAL distributions 
shown in Fig. Efb) and (c), which are also sensitive to the prompt-photon events. On 
average, the AE should be larger for the dijet events, where more energy is radiated 
outside the dijet system than for the 7+jet events. After the inclusion of the background, 
this distribution is well reproduced by MC. 

The background fraction described above was used in the calculation of the total 7+jet 
cross section. For differential cross sections, the background fractions were determined 
by fitting the BPRE signal independently in each bin of the respective distributions. In 
order to reduce statistical fluctuations in regions of small statistics, it was assumed that 
the background fractions varied smoothly from bin to bin. Therefore, the dependence of 
the background fraction on Ej,, vp ', E^ , rf ct and x° bs was obtained by fitting the back- 
ground fractions for each bin with a linear function. The number of prompt-photon events 
in each kinematic bin was determined from such a linear- regression fit. The statistical 
uncertainties on the number of signal events were evaluated using 68% confidence-level 
limits on the linear fit of the fractions. 

6 Cross section calculations and systematic 
uncertainties 

The differential cross sections for a given observable Y were determined as: 

da _ N 

dY ~ A-C-AY ' 

where N is the number of prompt-photon events in a bin of size AY, A is the acceptance 
and £ is the integrated luminosity. The acceptance was calculated using MC from the 
ratio of the number of reconstructed events after the selection cuts to the number of 
generated events. 

The systematic uncertainties were evaluated by changing the selection and the analysis 
procedure. The contribution of each cut variation to the total cross section is given in 
parentheses as a percentage of the total cross section: 

• the calorimeter energy scale was changed by ±3% (li 1 1 7 %); 

• the transverse momentum cut and the 77 range for the photon and hadron jet were 
lowered (raised) independently by one a of the resolution. The systematic uncertainty 
due to the transverse-energy cut for the photon was found to be (-3'%%)- The largest 



8 



systematical uncertainty due to the transverse-energy cut on the jet was The 
systematic uncertainty associated with the variations in the pseudorapidity was small 
(±0.8%); 

• the uncertainty in the quantity of inactive material in the MC was estimated by varying 
the amount of inactive material by ±5% of a X in the GEANT-based correction factors 

• the cone radius for the determination of the BPRE signal was changed by ±0.1 units 
(7.2 ?%)■ A larger cone size leads to a larger leakage of hadronic energy into the photon, 
which is not well simulated in MC; 

• variations of the cuts on tjjb, ^vertex and on total missing transverse energy (±2%); 

• the resolved contribution in MC was changed by ±15% (< 1%); 

• the cut on the electromagnetic fraction ^EMc/^tot f° r the photon jet was changed by 
±0.02 (l°j%); 

• the acceptance correction and the fraction of background photoproduction events was 
determined using Herwig (—0.6%). 

The overall systematic uncertainty was determined by adding the above uncertainties in 
quadrature. A 2% normalisation uncertainty due to the luminosity measurement error 
was not included in the systematic uncertainties. 

As an additional check, the differential 7+jet cross sections were found by using the global 
background fraction determined in the full kinematic range. Further, the cross sections 
were calculated from the number of the detector-level events in the data and MC after 
requiring a BPRE signal < 7 mip, i.e. in a region where the purity of the prompt photon 
sample is expected to be above 50%. The results from these alternative methods were 
consistent with the final cross sections. 

7 Results 

The total cross section for the process ep — > e + 7 pr ompt + jet + X for 0.2 < y < 0.8, 
Q 2 < 1 GeV 2 , 5 < E? T < 16 GeV, 6 < < 17 GeV, -0.74 < rf < 1.1, -1.6 < rf ct < 2.4 

was measured to be 

a(ep — > e + 7 P rom P t + jet + X) = 33.1 ± 3.0 (stat.) ^(syst.) pb. 

This cross section should be compared to the QCD predictions after the hadronisation 
corrections: 23.3+i;? pb (KZ), 23.5j££pb (FGH) and 30.7lfjpb (LZ). The scale uncer- 
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tainties on the QCD calculations were estimated by varying fi R between /xr/2 and 2/xr. 
The Pythia and Herwig cross sections are 20.0 pb and 13.5 pb, respectively. 

The differential cross sections as functions of and rj for the prompt-photon candidates 
and for the accompanying jets are shown in Figs. El and El Figure [7| shows the distribution 
for x° bs defined as J^vjetC^ — ^z) / (2E e y) (the sum runs over the photon candidate and 
the hadronic jet). Table ^ gives the differential cross sections with the statistical and 
systematical uncertainties, as well as the hadronisation-correction factors calculated in 
the same bins as the data 2 . 

The Pythia and Herwig differential cross sections do not rise as steeply at low Ej. 
as do the data. In addition, they underestimate the measured cross sections. The KZ 
NLO prediction, corrected for hadronisation effects as described in Sect. El describes the 
data better. However, it underestimates the observed cross section at low E^, and in 
the forward jet region. The observed difference between the data and the NLO QCD 
calculations is concentrated in the x° bs < 0.75 region which is sensitive to the resolved 
photon contribution. 

The FGH prediction is similar to the KZ NLO. The largest difference between the two 
predictions is found for the rf et cross section, where the FGH cross section is closer to the 
data in the forward jet region. The renormalisation scale uncertainty for the FGH QCD 
calculations is similar to that estimated for the KZ predictions (not shown). 

The LZ prediction based on the fey-factorisation approach corrected for hadronisation 
effects gives the best description of the Et and 77 cross sections. In particular, it describes 
the lowest Ej, region better than the KZ and FGH NLO predictions. The rf et cross section 
for the associated jet in the forward region is also better reproduced by the LZ calculation. 

It is difficult to compare the present cross sections with the HI result [8], since a sig- 
nificant model-dependent extrapolation to the low Ej? region used by HI is required. 
A comparison in the region E^ > 6 GeV shows good agreement with HI for the Ej? 
differential cross section. 

Since the largest difference between the NLO calculations and the data is concentrated in 
the region of low Ej, and low E^ , it is instructive to verify the level of agreement with 
NLO when the minimum transverse energy of the detected prompt photons is increased 
from 5 GeV to 7 GeV. In this case, hadronisation corrections are expected to be smaller. 
Further, in comparison with the previous measurements, such a choice may emphasize 
different aspect of contributions of high-order QCD radiation [36], since the transverse 
energy of the prompt-photon is larger than that of the jet. 

2 The actual hadronisation corrections applied to the NLO calculations shown in Figs. l5l9l were calculated 
using finer bins. 
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The total ep — > e + 7 pr0 mpt + jet + X cross section for Ej, > 7 GeV (keeping the other 
cuts the same as before) is a = 13.8 ± 1.2 (stat.) lj;'g(syst.) pb. This result agrees with 
the QCD calculations after the hadronisation corrections: 14.9^'gpb (KZ), 13.4*ggpb 
(FGH) and 13.6±?;g pb (LZ). The PYTHIA and HERWIG models predict 13.7 pb and 
9.4 pb, respectively. 

Figures |H1 El and Table |21 show the corresponding differential cross sections. The applied 
hadronisation corrections are given in Tabled For the Ej, > 7 GeV cut, both the NLO 
QCD and the LZ predictions agree well with the data. The Pythia MC model also agrees 
well with the cross sections, while Herwig is still below the data. 

8 Conclusions 

The photoproduction of prompt photons, together with an accompanying jet, has been 
measured in ep collisions at a centre-of-mass energy of 318 GeV with the ZEUS detector 
at HERA using an integrated luminosity of 77 pb -1 . 

In the kinematic region Ej, > 5 GeV and E^ > 6 GeV the prompt-photon data disagree 
with the available MC predictions which predict a less steep rise of the cross sections 
with decreasing Ej,. The discrepancy is reduced for the NLO calculations. However, they 
still underestimate the data in the low Ej, and Ej? regions, which are likely to be the 
most sensitive to the treatment of high-order QCD terms and hadronisation effects. The 
best description of the data was found for the calculations based on the /^-factorisation 
approach and unintegrated parton densities. 

When the minimum transverse energy of prompt photons is increased from 5 GeV to 
7 GeV, both NLO QCD and the /c^-factorisation calculations describe the data well. 

Acknowledgements 

We thank the DESY Directorate for their strong support and encouragement. The re- 
markable achievements of the HERA machine group were essential for the successful 
completion of this work and are greatly appreciated. We are grateful for the support 
of the DESY computing and network services. The design, construction and installa- 
tion of the ZEUS detector have been made possible owing to the ingenuity and effort of 
many people from DESY and home institutes who are not listed as authors. We thank 
M. Fontannaz, G. Heinrich, M. Krawczyk, A. Lipatov and A. Zembrzuski for discussions 
and for providing the QCD calculations. 



11 



References 



[1] L.E. Gordon and W.Vogelsang, Phys. Rev. D 52, 58 (1995). 

[2] M. Fontannaz, J. P. Guillet and G. Heinrich, Eur. Phys. J. C 21, 303 (2001); 
M. Fontannaz and G. Heinrich, Eur. Phys. J. C 34, 191 (2004). 

[3] M. Krawczyk and A. Zembrzuski, Phys. Rev. D 64, 14017 (2001); 

A. Zembrzuski and M. Krawczyk, Preprint hep-ph/0309308, 2003. 

[4] A.V. Lipatov and N.P Zotov, Phys. Rev. D 72, 054002 (2005). 

[5] ZEUS Coll., J. Breitweg et al, Phys. Lett. B 413, 201 (1997). 

[6] ZEUS Coll., S. Chekanov et al, Phys. Lett. B 511, 19 (2001). 

[7] ZEUS Coll., J. Breitweg et al., Phys. Lett. B 472, 175 (2000). 

[8] HI Coll., A. Aktas et al., Eur. Phys. J. C 38, 437 (2005). 

[9] ZEUS Coll., U. Holm (ed.), The ZEUS Detector. Status Report (unpublished), 
DESY (1993), available on http://www-zeus.desy.de/bluebook/bluebook.html. 

[10] N. Harnew et al., Nucl. Inst. Meth. A 279, 290 (1989); 

B. Foster et al, Nucl. Phys. Proc. Suppl. B 32, 181 (1993); 
B. Foster et al, Nucl. Inst. Meth. A 338, 254 (1994). 

[11] M. Derrick et al, Nucl. Inst. Meth. A 309, 77 (1991); 
A. Andresen et al., Nucl. Inst. Meth. A 309, 101 (1991); 
A. Caldwell et al, Nucl. Inst. Meth. A 321, 356 (1992); 
A. Bernstein et al, Nucl. Inst. Meth. A 336, 23 (1993). 

[12] J. Andruszkow et al., Preprint DESY-92-066, DESY, 1992; 
ZEUS Coll., M. Derrick et al., Z. Phys. C 63, 391 (1994); 
J. Andruszkow et al., Acta Phys. Pol. B 32, 2025 (2001). 

[13] S. Magill and S. Chekanov, Proceedings of the IX Int. Conference on Calorimetry 
(Annecy, Oct 9-14, 2000), B. Aubert et al. (ed.), p. 625. Frascati Physics Series 21, 
Annecy, France (2001). 

[14] V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972); 
L.N. Lipatov, Sov. J. Nucl. Phys. 20, 94 (1975); 
Yu.L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977); 
G. Altarelli and G. Parisi, Nucl. Phys. B 126, 298 (1977). 

[15] M. Gliick, E. Reya and A. Vogt, Z. Phys. C 67, 433 (1995). 

[16] M. Gliick, E. Reya and A. Vogt, Phys. Rev. D 45, 3986 (1992); 
M. Gliick, E. Reya and A. Vogt, Phys. Rev. D 46, 1973 (1992). 

12 



[17] M. Gliick, E. Reya and A. Vogt, Phys. Rev. D 48, 116 (1993). 

[18] A.D. Martin et al., Eur. Phys. J. C 14, 133 (2000); 
A.D. Martin et al, Eur. Phys. J. C 23, 73 (2002). 

[19] P. Aurenche, J.P Guillet and M. Fontannaz, Z. Phys. C 64, 621 (1994). 

[20] E.M. Levin et al, Sov. J. Nucl. Phys. 53, 657 (1991); 

S. Catani, M. Ciafaloni and F. Hautmann, Nucl. Phys. B 366, 135 (1991); 
J.C. Collins and R.K. Ellis, Nucl. Phys. B 360, 3 (1991). 

[21] M.A. Kimber, A.D. Martin and M.G. Ryskin, Phys. Rev. D63, 114027 (2001); 
G. Watt, A.D. Martin and M.G. Ryskin, Eur. Phys. J. C31, 73 (2003). 

[22] S.D. Ellis and D.E. Soper, Phys. Rev. D 48, 3160 (1993); 
S. Catani et al., Nucl. Phys. B 406, 187 (1993). 

[23] A. Zembrzuski. Private communication. 

[24] M. Bengtsson and T. Sjostrand, Comp. Phys. Comm. 46, 43 (1987). 

[25] T. Sjostrand, Comp. Phys. Comm. 82, 74 (1994). 

[26] T. Sj6strand et al., Comp. Phys. Comm. 135, 238 (2001). 

[27] G. Marchesini et al., Comp. Phys. Comm. 67, 465 (1992). 

[28] G. Corcella et al., JHEP 1, 10 (2001). 

[29] CTEQ Coll., H.L. Lai et al., Eur. Phys. J. C 12, 375 (2000). 

[30] G.A. Schuler and T. Sjostrand, Phys. Lett. B 376, 193 (1996). 

[31] R. Brun et al, geant3, Technical Report CERN-DD/EE/84-1, CERN, 1987. 

[32] F. Jacquet and A. Blondel, Proceedings of the Study for an ep Facility for Europe, 
U. Amaldi (ed.), p. 391. Hamburg, Germany (1979). Also in preprint DESY 79/48. 

[33] G.M. Briskin, Diffractive Dissociation in ep Deep Inelastic Scattering. Ph.D. 
Thesis, Tel Aviv University, 1998. 

[34] J.E. Huth et al., Research Directions for the Decade. Proceedings of Summer Study 
on High Energy Physics, 1990, E.L. Berger (ed.), p. 134. World Scientific (1992). 
Also in preprint FERMILAB-CONF-90-249-E. 

[35] S. Frixione and G. Ridolfi, Nucl. Phys. B 507, 315 (1997). 

[36] M. Fontannaz, J.P. Guillet and G. Heinrich, Eur. Phys. J. C 22, 303 (2001). 

[37] ZEUS Coll., S. Chekanov et al, Phys. Lett. B 573, 46 (2003). 

[38] P.R.B. Saull, A Monte Carlo Generator for Deeply Virtual Compton Scattering at 
HERA, 1999, available on http://www-zeus.desy.de/physics/diff/pub/MC. 



13 



(GeV) 


da/dEl (pb/GeV) 
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5.00,7.00 


9.0 il.llj;? 
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7.00,9.00 


3.7 ± 0.6i|j;jj 


1.01 


9.00,11.00 


2.9 ±0.51^ 


1.05 


11.00,13.00 


0.7±0.2lgj 


1.06 


13.00,16.00 
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1.06 


rp 


da/drf (pb) 




-0.74, -0.34 
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0.90 




da/drf et (pb) 




-1.60,-0.80 
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-0.80,0.00 


11.7 ± 1.4±J;£ 
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1.09 


™obs 
X 7 


da/dx° bs (pb) 




0.00,0.25 


3.9 ± 3.111a 


0.91 


0.25,0.50 


37.9 ± 7.9l? 2 8 3 


0.95 


0.50,0.75 


24.5 ± 5.21^ 


1.06 


0.75,1.00 


80.4 ± 7.2l? 2 4 6 


0.90 



Table 1: The differential prompt-photon cross sections with additional jet re- 
quirement measured in the region 0.2 < y < 0.8, Q 2 < 1 GeV 2 , 5 < E^ < 16 GeV, 
6 < E 3 T et < 17 GeV, -0.74 < rf < 1.1 and -1.6 < r] jet < 2.4. The statistical 
and systematical uncertainties are given separately. The hadronisation correction 
factors (see the text) applied to the QCD calculations for the same kinematic bins 
as for the data are also shown. 
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£fhad 
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0.25,0.50 


7.3 ± 1.51a;? 


1.11 
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9.1 ± 1.61^ 


1.16 
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33.9±4.4iy 


0.99 



Table 2: The differential prompt-photon cross sections with additional jet require- 
ment measured in the region defined as for TableU\ except for the cut on the trans- 
verse energy of the prompt photon, which was increased to 7 GeV . The statistical 
and systematical uncertainties are given separately. The hadronisation correction 
factors applied to the QCD calculations for the same kinematic bins as for the data 
are also shown. 
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Figure 1: Examples of diagrams for 7-/- jet events at leading order: direct photon 
interactions, (a-b) and resolved photon interactions, (c-d). The resolved diagrams 
with t-channel exchange are not shown. 
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Figure 2: The response of the BPRE detector to isolated photons in the DVCS 
data sample. The DVCS Monte Carlo distribution was normalised to the data. 
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Figure 3: The BPRE response to isolated photons, tt° and rj in the MC simulation. 
An initial energy of 10 GeV for all particles was used. The amount of inactive 
material in front of the BPRE was set to 1.25X . 
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Figure 4: Comparisons between the data and MC for: (a) the BPRE signal 
for the photon candidate; (b) the difference between the total calorimeter energy 
and the energies of the jet and the photon candidate, AE = E tot — E }et — E 1 ; 
and (c) the distance from the photon candidate to any EFO (see the text). The 
non-hatched histogram is the sum of the prompt photon MC and the background 
MC. The fraction of the prompt-photon events was found after a x 2 minimisation 
procedure for the BPRE distribution shown in (a). 
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Figure 5: The 7 + jet differential cross sections as functions of and 77 7 com- 
pared to theoretical QCD calculations (with hadronisation corrections included). 
The histograms show the predictions of the Monte Carlo models. The inner error 
bars show the statistical uncertainties, the outer ones show statistical and system- 
atic uncertainties added in quadrature. The shaded bands for the KZ prediction 
correspond to the uncertainty in the renormalisation scale which was changed by 
a factor of 0. 5 and 2. A similar uncertainty exists for the FGH prediction ( not 
shown). 
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Figure 6: The^+]et differential cross sections as functions of andrf ct com- 
pared to the QCD calculations (with hadronisation corrections) and Monte Carlo 
models. 
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Figure 7: The x° cross section for 7 + jet events compared to the NLO QCD 
calculations (with hadronisation corrections) and Monte Carlo models. 
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Figure 8: The differential 7 + jet cross sections as functions of: a) rp ', b) 
and c) rf ct compared to the QCD calculations (with hadronisation corrections) and 
Monte Carlo models. The cuts are the same as for Figs. [3] and 01 except for the 
cut on the transverse energy of the prompt photons, which was increased to 7 GeV. 
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Figure 9: The x° cross section for 7 + jet events compared to the QCD calcula- 
tions (with hadronisation corrections) and Monte Carlo models. The cuts are the 
same as for Fig. ^ except for the cut on the transverse energy of the prompt photon 
which was increased to 7 GeV . 
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